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Abstract: We formulate the string field theory in zero-dimensional target space 
corresponding to the two-dimensional quantum gravity theory defined through 
Causal Dynamical Triangulations. This third quantization of the quantum gravity 
theory allows us in principle to calculate the transition amplitudes of processes in 
which the topology of space changes in time, and to include non-trivial topologies 
of space-time. We formulate the corresponding Dyson-Schwinger equations and 
illustrate how they can be solved iteratively. 
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1. Introduction 

The formulation of a non-critical string theory in which the conformal mode 
plays an important role dates back to Polyakov. He emphasized the worldsheet 
formulation of string theory as a two-dimensional quantum gravity theory coupled 
to matter M. This triggered non-perturbative definitions of non-critical string 
theory 0, [3], M, introducing what is now called dynamical triangulation (DT) as 
a regularization of the worldsheet theory. When the dimension of space-time was 
larger than 1 these attempts in some sense did not work. One could show that 
the outcome was not a proper string theory, but a theory where the worldsheet 
had degenerated into branched polymers 0. However, when considering matter 
fields with central charge c < 1, these regularized theories led to what is now 
known as non-critical string theory, a very useful toy model of real string theory. 
In particular, it has been possible to formulate a string field theory of non-critical 
string theory |], [7[ which is very much simpler than the critical string field theory. 
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The use of causal dynamical triangulations (CDT) rather than DT 
ularization of quantum gravity was inspired by earlier ideas in one insists, 
starting from a Lorentzian space-time, that only causal histories contribute to the 
quantum gravitational path integral. In addition, one assumes the presence of 
a global time-foliation. In this way the space-times appearing in the regularized 
path integral become a set of piecewise linear causal geometries, made out of 
triangles (two-simplices) whose edge lengths provide an ultraviolet cut-off. For a 
detailed description of how to construct these geometries we refer to 0, |IU| in two 
dimensions and |TT] in higher dimensions. In order to perform the summation 
over these causal geometries we perform a rotation to Euclidean space-times. 
Each piecewise linear causal geometry as defined in [11 has a continuation to 
Euclidean signature, but the class of Euclidean geometries included in the path 
integral will only be a subclass of the total class of Euclidean geometries, and 
the result of the summation will therefore be different from that of Euclidean 
quantum gravity. 

One is interested in the limit where the lattice spacing a goes to zero. There 
is evidence for the existence of an underlying (non-perturbatively defined) con- 
tinuum quantum field theory in four dimensions |12|] and the results seem to be in 



qualitative agreement with recent renormalization group calculations [|13|] . These 



intriguing developments in the /bt/r-dimensional theory are based on numerical 
simulations, since analytical tools are presently unavailable. In two dimensions 
the situation is different, since the quantum gravity model can be solved analyt- 
ically at the discretized level and the limit a — > can be constructed. 

In we showed that the original two-dimensional CDT model of quantum 
gravity defined and solved in || can be generalized to a model where one al- 
lows for the creation of so-called baby universes, branching off from the "parent 
universe". The creation of a baby universe results in at least one point where 
from a Lorentzian point of view the metric is degenerate [15|. One cannot in- 
voke the classical theory to decide a priori whether or not such geometries should 
be included in the path integral. In |J we made the choice to suppress these 
configurations. We could also show that if they were completely unsuppressed 
one would recover Euclidean 2d quantum gravity as defined via DT or quantum 
Liouville theory. The converse was demonstrated in Jl6| : if one integrates out all 
baby universes in Euclidean quantum gravity, one obtains CDT. 

Quite surprisingly, there exists yet a third possibility, namely, a double-scaling 
limit where the creation of baby universes in CDT can be associated with the 
gravitational coupling constant flH |. In this double-scaling limit one can cal- 
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culate the disc amplitude and finds a result which is analytically connected to 
the old CDT result, the expansion parameter being the gravitational coupling 
constant. However, this cannot - at least not by simple analytic continuation - 
be connected to the Euclidean theory. Thus we have arrived at a theory which 
allows the creation of baby universes, but in a much more controlled way than 
in Euclidean quantum gravity. Of course, unlike the original CDT prescription, 
this construction contains causality-violating features at the level of the piecewise 
linear Lorentzian geometries. However, as we will see, the Lorentzian structure 
still plays a role in "taming" them. - Apart from the interesting observation 
that such a new theory exists, it may have important implications for the higher- 
dimensional theories. The attempt to formulate Euclidean higher-dimensional 
quantum gravity theories using DT as a regularization ran into the problem that 
baby universes completely dominate the path integral and make it difficult to 
obtain a physically sensible continuum limit. Now we see that there may exist 
a way to include the creation of baby universes in a controlled manner, starting 
with the CDT regularization of the quantum gravity theory. 

In this paper we show that the construction of [14] can be turned into a 



full-fledged third quantization of 2d quantum gravity. In the terminology of || 
this is a string field theory for c = 0, in the sense that it allows the calculation 
of amplitudes for splitting and joining of (spatial) universes and as well as the 
inclusion of different space-time topologies. 

The remainder of this article is organized as follows: In Sec. § we review 
briefly the results of the generalized CDT model. In Sec. |3] we show how to 
define a string field theory, and in Sec. |] we show how it reproduces the results 
of the generalized CDT model. In Sec. [5] we derive the general Dyson- Schwinger 
equations and in Sec. [5] we show how they can be used to calculate in a systematic 
way multi-universe and topology-changing amplitudes. Finally, we discuss the 
interpretation and possible generalizations in Sec. |7]. 

2. Generalized Causal Dynamical Triangulation in 2d 

We will initially assume that the two-dimensional space-time has topology S 1 x 
[0, 1]. After rotation to Euclidean signature, the pure gravity action is given by 

= A J d^det^C) + x£dl 1 + y£ d/ 2 , (2.1) 

where A is the cosmological constant, x and y are two so-called boundary cos- 
mological constants, is the metric of a geometry of the kind described above, 
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and the line integrals refer to the lengths of the in- and out-boundaries induced 
by g^y. The so-called •proper-time propagator is defined by 

G x (x,y-t) = J V[g^)e- S ^\ (2.2) 

This represents the Euclideanization of a functional integral over space-times with 
Lorentzian signature, performed over all causal geometries [g^] such that the 
final (or "exit") boundary with boundary cosmological constant y is separated 1 
a geodesic distance t from the initial (or "entrance") boundary with boundary 
cosmological constant x. To arrive at the integral ( |2.2[ ), all causal geometries 
have been rotated to Euclidean signature, a procedure which is well defined in 
the CDT regularization of the path integral. 

Calculating the path integral (|2.2j ) with the help of the CDT regularization 
and taking the continuum limit as the edge length a of the triangles goes to zero 
leads to the equation 

jG x {x, y; t) = - A [(^ _ X)Gx{Xi y . t) ] ; (2 .3) 

which can readily be solved ||. Let l\ denote the length of the initial and I2 
the length of the final boundary. Rather than considering a situation where the 
boundary cosmological constant x is fixed, we will take l\ as fixed, and denote the 
corresponding propagator by G\(li,y;t), with similar definitions for G\(x, fat) 
and G\(li, fat). All of them are related by Laplace transformations, for instance, 

/*oo /*oo 

G x (x,y;t) = dl 2 dhG x (h,fat)e- xh -y l \ (2.4) 



JO JO 

Next, we will turn our attention to the so-called disc amplitude, associated 
with a piece of space-time which has the topology of a disc. Strictly speaking, 
the disc amplitude does not exist in CDT. A spatial slice in a two-dimensional 
Lorentzian space-time of the type we are considering will by construction be a 
one-dimensional space-like subspace of topology S 1 , i.e. a circle. Now, there is 



1 The statement that the exit boundary is separated by a geodesic distance t from the entrance 
boundary means in this context that all points on the exit boundary have a geodesic distance 
t to the entrance boundary. The geodesic distance of a point on the exit loop to the entrance 
loop is defined as the minimal geodesic distance from the exit point to points on the entrance 
loop. In the piecewise flat, triangulated geometries we are working with, "distance" is given by 
"link distance" . In the case of the original, "pure" CDT without any causality violations, the 
notion of distance between boundaries just introduced is symmetric under exchange of entrance 
and exit boundary. For the generalized models discussed below, this will no longer be the case. 
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Figure 1: In all four graphs, the geodesic distance from the final to the initial loop 
is given by t. Differentiating with respect to t leads to eq. (|2.6|). Shaded parts of 
graphs represent the full, g-dependent propagator G\ ig and disc amplitude W\ %g , and 
non-shaded parts the CDT propagator G\. 



no way this can be extended to a well-defined Lorentzian geometry everywhere in 
the interior of any finite disc whose boundary is the circle. The light cones of the 
geometry must degenerate in at least a point, because the disc does not extend 
infinitely in time. However, after rotation to Euclidean signature 2 , we can define 
a disc amplitude, which is related to G\(x, l 2 ] t) by 

f°° 1 

W x {x)= dtG x (x,l 2 = 0;t) = =. (2.5) 

Jo 1 + vA 

There is clearly a latest time t where the spatial universe contracts to length zero 
and vanishes into the "vacuum" . When introducing the string field theory below, 
we will see that this process has a natural realization as a tadpole term in the 
string field Hamiltonian. 

We will now allow for the possibility that space branches into disconnected 
parts as a function of proper time t, and introduce a coupling constant g of mass 



2 We should emphasize that there is in principle a choice involved when generalizing the 
unique Wick rotation of CDT || to situations where the causal structure of the piecewise 

flat geometries has singularities. One might attach certain complex and/or singular weights to 
such singularities in the Euclideanization, for example, of the kind envisaged in [^5). When the 
disc amplitude was first introduced in two-dimensional CDT in order to compare it to Euclidean 
models g], no extra weight was associated with it, leading to the disc amplitude (|2.5|). In the 
present work, following [jf^, we will associate finite, real weights with baby universes and 
branching points, as will be explained in detail below. 
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Figure 2: Graphical illustration of eq. (|2.7|) . Shaded parts represent the full disc 
amplitude W\ <g , unshaded parts the CDT disc amplitude W\ and propagator G\. 



dimension 3 associated with the branching 3 . As shown in |L4|, this modifies the 
equation for the proper-time propagator to 



-G Ai5 (x, y- 1) = -— [((a; 2 -\)+2g Wy g (x) ) G x . (l (.r. //: I ) 



(2.6) 



where the generalized nature of the propagator G\ t9 is indicated by the additional 
subscript g. The graphical representation of the integral version of eq. ( |2.6| ) is 
shown in Fig. [!]. At this point, the new, generalized disc amplitude W\ >g (x) is 
unknown and has to satisfy the equation 



oo oo 

Wi%)+g J dt J &h&l 2 {h + l2)Gfl{xM + h]t)W x , g {h)W x>g {l 2 ) 



(2.7) 

where superscripts (0) indicate the CDT amplitudes introduced in eqs. (|2.5|) and 
( p.3|) above, that is, 



W%(x) = W Xa= o(x) = W x (x) 



and similarly for G^ . The graphical representation of eq. ( 2.7| ) is shown in Fig. 



|. The integrations in (2.7) can be performed and one finds |1 



3 One could in principle have considered a more general branching process, where more than 
one baby universe can sprout at any given time step t. However, starting with the discretized 
theory and a lattice cut off, one can show that such processes are suppressed when the lattice 



spacing goes to zero [14 . This is related to the fact that g has mass dimension 3. 
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where 



(x — c 



x + c) 2 — 



•2g 



u 



u' 



u + 



A3/2 



W x , 9 (x) = (x 2 -\) + 2gWx, g (x). 
Using the definition fl2.10|) and eq. ( [2.9D , we can write eq. 

d d r - 

-Q t G\, g {x, y; t) = -— \Wx t g{x) G x ,g(x, y\ t) 



as 



The solution is 



G x , g {x,y;t) 



W x , g {x{t,x)) 



W x ,g(x) x(t,x)+y 
where x) is the solution of the characteristic equation for ( |2.11| ), 



dx ~ 

-^ = -W X , 9 (x), X(0,X) 



X. 



(2.9) 
(2.10) 

(2.11) 
(2.12) 

(2.13) 



The generalized CDT model of 2d quantum gravity we have defined above 
is a perturbative deformation of the original model in the sense that it has a 
convergent power expansion of the form 



9 



.A 3 / 2 



(2.14) 



in the dimensionless coupling constant g/\ 3 ^ 2 . This implies in particular that the 
average number (n) of baby universes created during the proper-time evolution of 
the two-dimensional universe described by this model is finite, a property already 
observed in previous 2d models with topology change [ 17| . The expectation value 
of the number n of branchings can be computed according to 



<n> 



g dW x ,g(x) 
Wx, g (x) dg 



;2.i5) 



which is finite as long as we are in the range of convergence of W\^ g {x). This coin- 
cides precisely with the range where the function W\ )B (x) behaves in a physically 
acceptable way, namely, Wx, g {l) goes to zero if the length I of the boundary loop 
goes to infinity 



Is it possible to give a gravitational interpretation of the new coupling con- 
stant gl From a purely Euclidean point of view all graphs appearing in Fig. |2| 
have the fixed topology of a disc. However, from a Lorentzian point of view, which 
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comes with a notion of time, it is clear that the branching of a baby universe is 
associated with a change of the spatial topology, a singular process in a Lorentzian 



space-time ||15| . One way of keeping track of this in a Wick-rotated, Euclidean 
picture is as follows. Since each time a baby universe branches off it also has to 
end somewhere, we may think of marking the resulting "tip" with a puncture. 
(Of course, these baby universes can in turn have baby universes branching off 
them, giving rise to additional branchings and punctures.) From a gravitational 
viewpoint, each new puncture corresponds to a topology change and receives a 
weight 1/Gn, where Gn is Newton's constant, because it will lead to a change by 
precisely this amount in the two-dimensional (Euclidean) Einstein- Hilbert action 

Seh = -r4r~ / ^V9R- (2-16) 
2nG N J 

Let us view the continuum theory as the limit of a lattice theory (CDT) with 
lattice spacing a. On the lattice we have a dimensionless "bare" coupling constant 
<7o( Q ) — 9 a3 , where a is the lattice spacing (see |T^] for a detailed discussion). 
According to the arguments above we can now make the identification go(a) = 
g-i/Givw^ where Gn(cl) denotes the "bare" gravitational coupling constant. One 
can introduce a renormalized gravitational coupling constant by 

-J— = — 4- + -lnAa 2 . (2.17) 

GT G N (a) 2 1 ; 

This implies that the bare gravitational coupling constant Gjv(o) goes to zero 
like 1/| In a 3 1 when the cut-off vanishes, a — > 0, in such a way that the product 
e 1 / ^™ fX 3 ' 2 is independent of the cut-off a. We can now identify 

e-i/cr = g/x s/2 (218) 

as the genuine coupling parameter in which we expand. 

This renormalization of the gravitational (or string) coupling constant is rem- 
iniscent of the famous double-scaling limit in non-critical string theory 4 ||18|| . In 
that case one also has g oc e^ 1 ^^" , the only difference being that relation ( [2.1 7|) 
is changed to 

— — = - l , . + -InAa 2 , (2.19) 

whence the partition function of non-critical string theory appears precisely as a 
function of the dimensionless coupling constant g/X 5 ^. 

4 It is called the double-scaling limit since from the point of view of the discretized theory 
it involves a simultaneous renormalization of the cosmological constant A and the gravitational 
coupling constant Gn- In this article we have already performed the renormalization of the 
cosmological constant. For details on this in the context of CDT we refer to H. 
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3. String Field Theory 



In quantum field theory, particles can be created and annihilated if the process 
does not violate any conservation laws of the theory. In string field theory, one 
has analogous operators which can create and annihilate strings. From the 2d 
quantum gravity point of view we are dealing with a third quantization of gravity: 
one-dimensional universes can be created and destroyed. In |J such a formalism 
was developed for non-critical strings in a zero-dimensional target space (or 2d 
Euclidean quantum gravity). We will follow the procedure outlined there closely 
and develop a string field theory or third quantization for CDT, which will allow 
us in principle to calculate any amplitude involving the creation or annihilation 
of universes. 

As starting point we assume the existence of a vacuum from which universes 
can be created. We denote this state |0) and define creation and annihilation 
operators through 

[*(0, * f (0] = «(* - 0. *(0I«> = (o|* f (0 = °- (3- 1 ) 

This assignment corresponds to working with spatial universes where a point has 
been marked. This is merely a formal aspect, to avoid having to put in certain 
combinatorial factors by hand when gluing universes together. The operators 
iff (I) and iff* (I) will be assigned the dimensions dim [iff] = dim [iff*] = 0. 

We could alternatively have chosen creation and annihilation operators which 
create and annihilate universes without such a mark. Instead of ( |3.1|) we then 
would have had 

mi), * f (oi = r^i - o, *(oio> = <oi*+(o = o, (3.2) 

with corresponding dimensional assignments dim [iff] = 1 and dim [iff * ] = 1. One 
could even let iff* create marked universes and iff annihilate unmarked universes if 
one compensated for the missing combinatorial factors by hand. In the following 
we will use the assignment ( ^.1[ ). 

Let us write the propagator equation (|2.3|) using the boundary length rather 
than the boundary cosmological constant as a variable 5 , that is, 

-G A (/ 1; l 2 ; t) =h(gp- X)G x {h, k] t), (3.3) 



5l 



For convenience of notation we have in (3.3) also marked the exit loop li in order to have 



symmetry between the loops at initial and final time, i.e. G\(h,l2',t) — l2G\(li,l2',t). 
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which we can also write as 

G x (li, h; t) = (h\e- tm) \h), H (l) = -l±L + xi (3.4) 

Associated with the spatial universe we have a Hilbert space on the positive 
half-line, and a corresponding scalar product 

<iMv*> = / jm)Mi), (3.5) 

which makes H (l) hermitian. The introduction of the operators \l/(7) and ^(l) 
in (|3.1|) can be thought of as analogous to the standard second quantization 
in many-body theory. The single-particle Hamiltonian becomes in our case the 
"single-universe" Hamiltonian H Q (l). It has normalized eigenfunctions ip n (l) with 
corresponding eigenvalues e n = 2n\f\ 1 n = 1,2,..., 

iP n {l) = Z e -^W_ X (Z), H (l)1> n (l) = e„V n (Z), (3.6) 

where p n -i(l) is a polynomial of order n — 1. We now introduce creation and 
annihilation operators a) n and a n corresponding to these states, acting on the 
Fock vacuum |0) and satisfying [a n , ajj = 5 nm . We define 

n n 

and from the orthonormality of the eigenfunctions with respect to the measure 
dl/ 1 we recover (|3.1|) . The "second-quantized" Hamiltonian is given by 

H = J o j^(l)H (l)*(l), (3.8) 

and the propagator G\(h, 1% t) is now obtained as 

G x (h,k;t) = (0|*(Z 2 )e-^ ^(Z 1 )|0). (3.9) 

While all of this is rather straightforward, the advantage of the formalism is that 
it automatically takes care of symmetry factors (like in the many-body applica- 
tions in statistical field theory), both when many spatial universes are at play 
and when they are joining and splitting. Following [H], we define the Hamiltonian 

H = H - gjdhj dl 2 ¥(h)^(l 2 )^(h + h) (3.10) 
-ag J dh J dh*\h + h)Hk)Hh)- J f P(0*(0> 
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Figure 3: The elementary terms of the string field theory Hamiltonian ( 3.10 ): (a) 
the single spatial universe propagator, (b) the term corresponding to splitting into two 
spatial universes, (c) the term corresponding to the merging of two spatial universes 
and (d) the tadpole term. 



describing the interaction between spatial universes (the different terms are illus- 
trated in Fig. |H). Here g is the coupling constant of mass dimension 3 we have 
already encountered in Sec. |||, and p(l) denotes the amplitude for a universe com- 
ponent of length I to disappear into the vacuum. The factor a has merely been 
introduced to distinguish between the action of the two terms proportional to g 
in ( p.lQ ) when expanding in powers of g. We will usually assume a = 1, unless 
explicitly stated otherwise. Note that the signs of all the interaction terms in 
( p.lO| ) are negative. This reflects the fact that we want these terms to represent 
the insertion of new geometric structures compared to the "free" propagation 
generated by H . These structures should therefore appear with positive weight 
when we expand e~ tH . 

The Hamiltonian H is hermitian except for the presence of the tadpole term. 
It tells us that universes can vanish, but not be created from nothing. The 
meaning of the two interaction terms is as follows. The first term replaces a 
universe of length li + 1 2 with two universes of length li and l 2 - This is precisely 
the process shown in Fig. [I]. The second term represents the opposite process 
where two spatial universes merge into one, i.e. the time-reversed picture. The 
coupling constant g is seen to appear as a string coupling constant: one factor of 
g for the splitting, and another factor of g for the merging of spatial universes, 
leading to a combined factor g 2 whenever a handle is added to the space-time. 

In a way the appearance of a tadpole term is more natural in CDT than in the 
original Euclidean framework in H], since in CDT it has its origin in a physical, 
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causality-violating process located at the end point (in time) of the disc, where 
the baby universe disappears into nothing, as we saw in Sec. |[ The tadpole term 
is a formal realization of this. Because of the one-to-one correspondence between 
punctures and baby universe branchings, we can also associate this process with 
the gravitational coupling constant, in this way linking it to g. The shift in 
associating the coupling g from the splitting of spatial universes to the vanishing 
of universes can be made explicit in our string field Hamiltonian H in ( |3.10| ). 
In ( |3.10p , the coupling constant g was associated with the splitting and joining 
of spatial universes, but no coupling constant with the tadpole term, i.e. the 
vanishing of a spatial universe. However, by redefining ^ and to 

# = & = g&, (3.11) 

the coupling constant g is shifted from the splitting to the tadpole term, i.e. 
precisely the shift mentioned above. In addition, the term associated with the 
joining of spatial universes will have the coupling constant g 2 , which precisely 
accounts for the change in topology. 

Finally, let us identify the true, dimensionless coupling constant governing 
( P-1C| ) . This can be done by re-expressing everything in units of 1 / y/X, which 
represents the natural length scale of our universe. Introducing the dimension- 
less length variable I = ly/X, the dimensionless boundary cosmological constant 
x = x/vA, the dimensionless time variable t = ty/X, the dimensionless tadpole 
density p(l) = p{l)/y/X, the dimensionless coupling constant g = g/X 3 ^ 2 (already 
introduced in eq. ( |2.14Q ), and finally the dimensionless Hamiltonian H = H/y/X, 
we can write 

ytft(0#o(0*(0, (3.12) 

where = \V(l) and \V^(l) = satisfy the same commutation relation as 

\V^(l) when expressed in terms of /, and H = V~XH, where 

H = H Q -g j dhj dW(h)&(h)Kh + ^) (3.13) 

-a~g J dh J dh^ih + l^ih^ih)- J |p(ZW). 

From this expression we can read off that the true coupling constant of the theory 
is the dimensionless quantity g, precisely the "double-scaling" coupling constant 
which already appeared in the calculation of W\ >g {x) and G\ ;g (x,y), c.f. eq. 
( |2.14| ). From the discussion above we also observe that the parameter associated 
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with a topological expansion of space-time is given by g 2 = g 2 /X 3 - In principle 
we can now calculate any process which starts from m spatial universes at time 
and ends with n universes at time t, represented by the amplitude 

G x , g (h, .., l m , l[, .., C t) = (Om'i) ■ ■ ■ *(0 e~^v]>t( Zl ) . . . *t( Zm )| ). (3.14) 

4. The a = limit 
4.1 The disc amplitude 

Let us now consider the simplest such amplitude, that of a single spatial universe 
disappearing into the vacuum. This is precisely the disc amplitude of generalized 
CDT considered in Sec. @. There, we allowed baby universes to branch off, but 
they were forbidden to rejoin the parent universe, and thus were destined to 
disappear into the vacuum eventually. In other words, the topology of space-time 
was not allowed to change during evolution. This can be incorporated in the 
string field-theoretic picture by choosing a = in ( p.lCj ). The disc amplitude can 
then be expressed as 

W X M) = lim W x M,t) = lim(0|e-^ (Q=0) ^ t (/)|0). (4.1) 

t— »oo ' t^oo 

It describes all possible ways in which a spatial loop can develop in time and 
disappear into the vacuum without changing the topology of space-time. Note 
that the tadpole term in ( |3.1(J| ) is needed if the amplitude ( |4.1| ) should be different 
from zero, since the state \l) = ip^(l)\0) is orthogonal to the vacuum state |0). 
We note that for a = the vacuum expectation value 

(0| e -^°)*t(i 1 )...*t( Zm )|o) = (4.2) 

(0| e-^( a=o) * t (Z 1 )|0)(0| e-^ (a=0) * t (Z 2 )|0)(0| • ■ • |0)(0| e~ tfl< - a=0) ^(l m )\0} 

factorizes, as one can easily prove using the algebra of the ^'s. This is an ex- 
pression of the fact that if we start out with m spatial universes, there is no way 
they can merge at any time if a = 0. 

Following ||, we obtain an equation for W\ >g (l) by differentiating fl4.1p with 
respect to t and using that H\0) = 0, 

= lim ^-W x , g (l,t) = lim(0|e-^ a=0 >[£(a = 0),¥(l)]\0). (4.3) 

The commutator can readily be calculated and after a Laplace transformation 
eq. (|4.3|) reads 

£ ((x 2 - X)W Xi9 (x) + gWl g {x)) = p(x), (4.4) 
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where the last term on the left-hand side of eq. (j4.4|) is a consequence of the fac- 
torization ( |4~2l) . Eq. ( |4.4|) has the generalized CDT solution ( |2 . 9|) - ( p7T0|) discussed 
in Sec. || if 

p(x) = l, i.e. p(l)=5(l). (4.5) 

This is a reasonable physical requirement, which we will implement in what fol- 
lows: the spatial universe can only vanish into the vacuum when the length of 
the universe goes to zero. 

4.2 Inclusive amplitudes 

After reproducing the generalized CDT disc amplitude W\ t9 (x) as the connected 
amplitude arising in the string field theory in the limit a = 0, we now want 
to understand how to rederive the proper-time propagator G\ g (x,y,t) in this 
context. This propagator is characterized by an entrance loop at time t = and 
an exit loop at time t, and also contains baby universes which branch off and can 
extend in time beyond time t, if only they vanish into the vacuum eventually, as 
indicated in Fig. [j]. 

We can reproduce this result in the a = limit of the string field theory 
by introducing the so-called "inclusive" Hamiltonian ||. Since we are working 
in the a = limit, universes can only branch and not merge during the time 
evolution, and all but one have to vanish into the vacuum. The branching process 
is associated with the term 

9 j dhj dl 2 ¥(h)¥(l 2 Mh + l 2 ) (4.6) 



in the Hamiltonian H of eq. ( |3.10| ). Once the branching has occurred, only one 
of the two universes can connect to the exit loop at time t, the other one has to 
continue until it eventually vanishes into the vacuum, a process which may occur 
at a time later than t. This scenario is captured by replacing 

&(h)^(l 2 ) - W x , 9 (h)^(h) + &(h)W x>g (l 2 ) (4.7) 
in eq. ( |4.6| ), thus arriving at the "inclusive Hamiltonian" 

Hind = j j^(l)H {lMl) - 29 J dh J dl 2 W x>g {h)^(l 2 Mh + l 2 ), (4.8) 

which enables us to rewrite the corresponding propagator G\ t9 (li, l 2 ; t) as 

G Xi9 (h,l 2 ;t) = (0|^(/ 2 )e-^-^ t (/ 1 )|0). (4.9) 
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Differentiating eq. (|4.9|) with respect to t, commuting Hi nc i through to the right 
and using H inc i\0) = 0, one obtains after a Laplace transformation eq. (|2^). 
We conclude that also the generalized CDT proper-time propagator has a simple 
string field-theoretic description. 

4.3 Propagator identities 

Our starting point was the functional integral (|2.2| ) over all two-dimensional ge- 
ometries with cylindrical topology, where the entrance and exit loop were sepa- 
rated by a geodesic distance t. This proper-time propagator played an important 
role, motivating the introduction of the string field Hamiltonian. As explained 
in footnote [I], this construction is based on a particular definition of the geodesic 
distance between the exit and the entrance loop: every point on the exit loop has 
geodesic distance t to the entrance loop, i.e. the minimal distance of a given point 
on the exit loop to the points on the entrance loop is precisely t, independent 
of the point on the exit loop. This implies that the exit loop as a whole has a 
specific, "parallel" orientation relative to the entrance loop. This is a very useful 
property for the propagator to have, ensuring the existence of simple composition 
rules and thus a Hamiltonian. 

What we will show next is that one can define more general amplitudes, 
which depend on a somewhat looser notion of distance between their boundary 
components, and which are obtained by appropriately gluing together proper- 
time propagators and disc amplitudes. These "combined" propagators obey non- 
trivial identities analogous to identities first found and verified in the Euclidean 
framework of the non-critical string field theory of ||. The fact that our CDT 
geometries still carry some memory of their original Lorentzian structure after 
mapping them to the Euclidean sector makes the physical interpretation of these 
identities in the CDT string field theory less clear, since the nature of the identities 
is rather "Euclidean" , as we shall see. 

The geometric configurations we are interested in consist of two entrance 
loops from which two universes propagate to the future, and then join to form 
a single universe, which eventually disappears into the vacuum. Two distinct 
configurations of this type are illustrated in Fig. [|. They differ in how much time 
elapses in each of the "legs" before they join. When summing over all geometries 
of fixed leg lengths (£1,^2), the legs will correspond to proper-time propagators 
of length t\ and £2, and the remainder of the geometry will correspond to a disc 
amplitude with boundary length I + 1', which has been pinched in a point such 
that it can be glued to the two exit loops of the propagators, of length I and 
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Figure 4: Two universes (whose time extensions add up to t) merging into a single 
one, and subsequently vanishing into the vacuum. The figure on the left shows the 
degenerate case where one leg has length t and the other length 0, whereas the figure 
on the right has two legs of unequal, non-zero length. An explicit computation shows 
that summing over all space-times of the first type gives the same result as summing 
over all space-times of the second type, for any choice of t\. 



V respectively. We will be interested in comparing situations where the two leg 
lengths sum to the same number t, such that = t — t\, for different t\. The 
left illustration in Fig. |4] corresponds to the extreme case t\ = 0, and the right 
one to some intermediate choice t\ < t%. We are not primarily concerned with 
the physicality or otherwise of these geometries, but simply note that they are 
well defined in our string field-theoretic set-up after Euclideanization, and possess 
calculable amplitudes. We can allow the propagation to be of the most general 
a = kind 6 , such that the dynamics is described by the inclusive Hamiltonian 
Hind- The two situations depicted in Fig. [| correspond to the two calculations 

POO 

/ dlG x , g (hJ;t)W x , g (l + l 2 ) (4.10) 
Jo 

POO POO 

/ dl / dl' G x>g {h, 1; h)W x , g {l + l')G Xjl (l 2 , 1'; t-h), (4.11) 
Jo Jo 



and 

f oo 



6 Strictly speaking, the processes described by ( 4.10 ) and ( 4.11 ) are of order a in the string 
field-theoretic terminology we have introduced above, since they describe two merging spatial 
universes. A related amplitude ^0(^1,^2), obtained by integrating over all times t, will be 
introduced later, c.f. (6.21). 
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for some < t\ < t. The remarkable fact is that the results of both calculations 
coincide! Equivalently, one can show that 



o rco l-OO 

° = ^T dl dl' G Xtg {hMi)W x ,S + l')Gx, 9 {k^^-ti)- (4-12) 
at i Jo Jo 

After a Laplace transformation, eq. ( 4.12| ) reads 

= 75- / —G^-z-hWx^Gx^-z-t-^). (4.13) 

Using the explicit form of G\, g (x, y; t), eq. ( |2.12| ), we can perform the ^-integration 
in eq. (|4.13| ). Next, with the help of eq. ( |2.10| ), we can express W x , g {x) in terms 
of W\, g (x) (given by eq. fl2.7|)), and finally, using eq. ( |2.13| ), we can perform the 
ti-differentiation. The result is 







o 2 



dxdy 



x\t ll x)-f{T-t ll y)) 1 (4.14) 



which is satisfied. The upshot of this calculation is that we can define a more 
general amplitude 

poo poo 

Q(h, l 2 , t) := / dl dl' G x , g {h, I; h)W Xt9 {l + l')G Xt9 {hi I'; t - h) (4.15) 
Jo Jo 

associated with this merger process, which only depends on the combined distance 
t along the legs, and which - as we have just proved - is invariant under how t is 
split into two. It is somewhat surprising that this invariance property holds, since 
as Lorentzian geometries the two situations depicted in Fig. f|are clearly distinct. 
Of course, during the Wick rotation the special character of the causality- violating 
merger point between the two "trouser legs" disappears, which may explain the 
validity of ( 4.10| ), just like in the Euclidean formulation 7 . In line with the latter, 



7 The corresponding equation in the case of non-critical string theory is 

d 2 







dxdy 

leading again to the result that the amplitude Q(li,l2',t) is independent of the subdivision of 
t = t\ + t-z- Note that in this purely Euclidean formulation a relation like ( 4.10 ) appears as 
a (necessary) consistency condition, whereas in the CDT case it is satisfied as a non-trivial 
identity. Moreover, we have in the Euclidean string field theory setting the additional consis- 
tency test that J* °° dtQ(h,h;t) = Q(h,fo), where Q(h,fo) is the so-called universal loop-loop 
correlator calculated from matrix models [[19l E0| . This was verified in 101 . 
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one may interpret the quantity G(h,h,t) as a generalized amplitude with two 
boundaries separated by a distance t, where the "distance" between two spatial 
loops is now defined as the smallest geodesic distance between any pair of points 
on the two loops, with no further constraints on the relative position of the two 
loops. In particular, this makes G{h,h,t) symmetric under the exchange of l\ 
and 1%. 



5. Dyson- Schwinger equations 

The disc amplitude W\ >g is one of a set of functions for which it is possible to 
derive Dyson-Schwinger equations. Here we will consider a more general class of 
functions. Defining the generating function Z(J;t) by 

Z(J;t) = (0\e- tfl ef dumHl) \0), (5.1) 

we have 



|e-'^t (i0 ...* t(( „ )|0)= _Z^«|_ 



(5.2) 

J=0 



For the special case of vanishing coupling a = 0, we have already seen that the 
amplitudes factorize, such that 

Z(J,t;a = 0) = e S diJ{i)w x Ji,t) ^ ( 53 ) 

where W\ >g (l,t) denotes the disc amplitude where the universe decays into the 
vacuum before or at time t, and where W\ >g (l,t = oo) is the disc amplitude we 
have already calculated. 

Following |J, we can obtain the Dyson-Schwinger equations in the same way 
as for the disc amplitude, the only difference being that when the constant a is 
no longer zero, these equations do not close but connect various amplitudes of 
more complicated topology. However, as we shall see, the equations can still be 
solved iteratively. We denote 

Z(J) = lim Z(J;t), (5.4) 

t— »oo 

Z(J) being the generating functional for universes disappearing into the vacuum. 
We now have 

= _ lim I- ( |e-^ e I^mm)\ ) = h m (0\e- tA Hef du ^^\0). (5.5) 

t^oo (jt t— >oo 
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Commuting the \I/(Z)'s in H past the source term effectively replaces these oper- 
ators by U(l), after which they can be moved to the left of any ^(l) and outside 
(0|. After that the remaining ^'(^'s in H can be replaced by 5/5J(l) and also 
moved outside (0|, leaving us with a integro-differential operator acting on Z(J), 

= J™ dlJ(l)o(l,J,-^)Z(J), (5.6) 



where 



°^ J -lj) = H ^lW)- m (57) 

s 2 , r ,„„ 5 



rl ^2 roc 

- 9l l dl, 5J(l')5J(l-V)- a9l l dl ' l ' J{1 ' 



SJ(l + V) 



The generating functional Z(J,t) also includes totally disconnected universes 
which never "interact" with each other. Since our main interest is in universes 
whose space-time is connected, the appropriate generating functional F(J,t) for 
connected universes is obtained by taking the logarithm of Z(J, t), following stan- 
dard quantum field-theoretic methods, 

F(J,t) = log Z(J,t). (5.8) 

From this we obtain the correlators 



(5.9) 



J=0 



and it is straightforward to translate the Dyson-Schwinger equation (|5.6|)-(|5T7D 
into an equation for the connected functional 

F(J) = lim F(J,t), (5.10) 

t— >oo 

namely, 

5 2 F(J) 



= l°°dl J(l) j#o(0 - - gl J dl> 



8J(l')8J(l - V) 



JF{J) 5F{J) r^rn* SF (J) 



-° l j/ l W)m^)- a9l h OT wr%r (5 - n) 

From eq. ( |5.11 ) one obtains the Dyson-Schwinger equation by differentiating 
( p. 11| ) with respect to J (I) a number of times and then taking J(l) = 0. 
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6. Application of the Dyson- Schwinger equation 

Let us introduce the notation 

5 n F(J) 



w(h, ...,l n ) 



5J(h)---5J(Q 
as well as the Laplace transform 



(6.1) 



j=o 



POO POO 

w{x x ,...,x n ) = dh--- dl n e- xlh -- Xnln w(h,...,l n ). (6.2) 
Jo Jo 



Next, differentiate eq. ( |5.11 ) with respect to J(l) one, two and three times, then 



take J(l) = and Laplace-transform the resulting equations. This leads to the 
following three equations (recall that H (x)f(x) = d x [(x 2 — X)f(x)}): 

= Hq(x)w{x) — 1 + gd x (w{x, x) + w(x)w(x)^ , (6.3) 

= (H (x) + H (y))w(x, y) + gd x w(x, x, y) + gd y w(x, y, y) (6.4) 

"w(x)— w(y)~ 



-2g (d x [w(x)w(x, y)]+d y [w(y)w(x } y)]) + 2agd x d y 



x-y 



= (H {x) + H (y) + H (z))w(x, y, z) (6.5) 
+gd x w(x, x, y, z) + gd y w(x, y, y, z) + gd z w(x, y, z, z) 
+2gd x [w(x)w(x, y, z)] + 2gd y [w(y)w(x, y, z)\ + 2gd z [w(z)w(x, y, z)) 
+2gd x [w(x, y)w(x, z)\ + 2gd y [w(x, y)w(y, z)\ + 2gd z [w(x, z)w(y, z)\ 

+2ag ( d x d y w{x ' Z) ~ w{y > Z) + d x d z W ^ V) ~ W ^ Z) + d y d z W ^ V) Z) 
\ x — y x — z y — z 

The general structure of these equations should now be clear 8 . We can solve the 
Dyson- Schwinger equations iteratively. To this end, introduce an expansion of 
w(xi, . . . , x n ) in powers of the coupling constants g and a, 

oo oo 

w(x 1 , . . . , x n ) = ^2 ak X] 9 m w^!, . . . ,x n ;m,k). (6.6) 

k=n— 1 m=k— 1 

The amplitude w(xi, . . . ,x n ) starts with the power (ctg)"™ 1 since we have to 
perform (n— 1) mergings during the time evolution in order to create a connected 



8 Interestingly, one can find a matrix model which reproduces the Dyson-Schwinger equations, 

m 
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geometry if we begin with n separated spatial loops. One can find the lowest-order 
contribution to w{x\) from (|6.3|) , use that to find the lowest-order contribution to 
w(xi, x 2 ) from (|6.4j ), and then use this again in (|6.5|) , which involves w(x\, x 2 , x 3 ), 
etc. Returning to eq. ( |6.3|) , we can use the lowest-order expression for w(xi, x 2 ) to 
find the next-order correction to w(x\), use this and the lowest-order correction 
for w(xi, x 2 , X3) to find the next-order correction to w(xi,x 2 ), etc. 

Two remarks are in order: first, the integration constants arising during 
the integration of ( |6.3|) -( |6~5l) and the corresponding higher-order equations are 
uniquely fixed by the requirement that the correlation functions fall off as the 
lengths li — > 00, i.e. the requirement that the Laplace-transformed amplitude 
w(xi, . . . ,x n ) be analytic for Xi > 0. Second, the expressions obtained for 
w(xi, . . . ,x n ) can of course be obtained directly from a diagrammatic expan- 
sion, using the interaction rules shown in Fig. ^, where the propagation is defined 
by Hq, and then integrating in a suitable way over the times tj involved. The 
results for the first few orders are 

w(x;0,0) = —^-=, (6.7) 
x + v A 

w(x;l,0)= x + 3v ^ , (6.8) 

w(x, y; 1, 1) = — -= J: , (6.9) 

1 ^ 2v / A(x + v / A) 2 (y + v / A) 2 

7A5 + 5A(x + y + z) + 3V\(xy + xz + yz) + xyz 
w(x,y,z;2,2) = = — — — — . 6.10) 

4A§ (VX + x) 3 (V\ + y) 3 (V\ + zf 

For all of these amplitudes, the space-time topology is trivial. To lowest order in 
g, i.e. without any additional baby universes, and using the results ( B.7|) -( |B~T0" 



111 



the iteration as described above, the genus-one and genus-two amplitudes become 



15A§ + HAx + 5y/\x 2 + 



?> 



w(x; 2, 1) = , (6.11) 

32A2(VA + a;) 5 

w(x;3,2) = j^—- fll319Ai + 19951A 3 a; + 21555Afa; 2 + 

2048A^(v / A+x) 9 V 

16955AV + 9765Afx 4 + 3885Ax 5 + 945v / Ax 6 + 105x 7N ) . (6.12) 



In a diagrammatic notation, the genus-two amplitude w(x;3,2) corresponds to 
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the following three diagrams (including suitable integrations over the times £j): 




w(x; 3, 2) 



As mentioned above, the expansion of the amplitude w(xi, . . . , x n ) starts with 
the power (ag) n_1 , coming from merging the n disconnected spatial universes. 
The remaining powers of a are associated with a non-trivial space-time topology 
in the form of h additional "handles" on the connected world sheet. From a purely 
Euclidean point of view this suggest a reorganization of the series according to 



W(Xi 



, ... j x n 



(ag) n l ^2(ag 2 ) h w h (x h ...,x n ), 



(6.13) 



h=0 



W h (Xl, 



y^^g j w(xi, . . . ,x n ;n — 1 + 2h + j,n— 1 +h), (6.14) 

j=0 



amounting to a topological expansion in ag 2 , solving at each order for all possible 
baby-universe creations which at some point will vanish into the vacuum. This 
implies that Wh(x\, . . . , x n ) will be a function of g, although we do not write the 
dependence explicitly. 

The Dyson- Schwinger equations allow us to obtain the topological expansion 
iteratively, in much the same way as in our earlier power expansion in g. Since 
we have w(x,x) = 0(a), this term does not contribute to lowest order; from eq. 
( |6.3D we obtain a closed equation for w (x), namely, 

H Q (x)w (x) + gd x wl(x) = I. (6.15) 

This equation is of course just eq. (fOj), where we have made the identification 

w (x) = W g>x (x). (6.16) 

Knowing w (x) allows us to obtain w (x, y) from (|6.4] ), since w(x, y, z) is of order 
0(a 2 ). Therefore the three- loop term does not contribute to the lowest order in 
a of eq. ( |6.4|) , which is 0(a), and we find that to lowest order 

w (x)-w (y) 
x-y 



H (x) + 2gd x w (x) + H (y) + 2gdyWo(y) )w (x,y) =-2d x d. ; 



(6.17) 
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We conclude that Wq(x, y) is entirely determined by the knowledge of wq{x). Note 
that using the definition fl2.10|) we can simplify ( |6.171 ) to 

— (Wx, g (x)w (x,y) )+— (WxMM^y)) = — Tnr 



(6.18) 

The solution w (x, y) can readily be found from eq. fl6.18| ), yielding 
I l 1 1 f {(x + c) + (y + e)} 2 \ 

7Mf¥)*<)\ 1/M + /(y)P "V- ( " 9) 

where 



f(x) = ^(x + c) 2 - 2g/c = Wx, g (x)/(x - c). (6.20) 
In fact, this solution was already found in WM since we have by definition that 



wo(x,y)= dtGx, g (x,y;t), (6.21) 

where Qx, g (x,y;t) is the Laplace transform of G(h,h]t) defined in ([4.111) , with 
t\ = t/2. This function is precisely the loop-loop function of |L4]]. When expanded 
to lowest order in g, it reproduces (|6.9|) . 

As should by now be clear, one can iterate the Dyson-Schwinger equations in 
a systematic way as a power series in the number h of handles of the world sheet, 
exactly like we iterated them as a function of the coupling constant g, leading to 

w(x) = wo(x) + ag 2 w\(x) + a 2 g W2(x) + • • • , (6.22) 
w(x, y) = agw (x, y) + a 2 g 3 w 1 (x, y) H , 

etc. As an instructive example we will calculate the genus-one amplitude w\(x). 
While eq. ( |6.15| ) was the 0th order in a of eq. (|6.3|), the 1st order reads 





y Wx, 9 (x) Wl (x) +w (x,x)) = 0, (6.23) 

where Wo(x, x) is given by eq. fl6.19| ). The integration constant is fixed by the 
requirement that w\{x) be analytic for x > 0, i.e. that w\(l) fall off as I —>■ oo. 
We obtain 

_ w (c, c) - w (x,x) _ (x + 3c) (x 2 + 2cx + 5c 2 - Ag/c) , , 

Wl[X) ~ W x , g (x) ~ 2c(4c 2 - 2^/c) 2 ((x + c) 2 - 2g/cfl 2 ' 1 ' j 

which upon expansion in powers of g to lowest order reproduces ( |6.11 ), as one 
would expect. 
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7. Discussion 



In the present work, we have developed a string field theory in zero-dimensional 
target space, based on the CDT quantization of two-dimensional quantum gravity. 
It shares many properties of the non-critical string field theory originally defined 
in ||, from which we borrowed the formalism in the first place. Yet, our results 
are different and in some ways simpler. The tadpole term in our case is simply 
p(l) = 5(1), encoding the fact that universes can only disappear into the vacuum 
if they have zero spatial volume (that is, zero length). This is in accordance with 
the interaction between spatial universes, which also preserves the total length. 
In non-critical string field theory the evolution in proper time results in a process 
where the original spatial universe at proper time t = spawns an infinity of 
(infinitesimal) baby universes during the time evolution. This is related to the 
fact that the proper time in non-critical string field theory has the anomalous 
length dimension 1/2. In our new CDT-based string field theory the situation is 
different. The proper time t has canonical dimension 1, and the number of baby 



universes created during the time evolution is finite [[LJ . 

It is not possible to connect the non-critical string field theory and the CDT- 
based string field theory by a simple analytic continuation in the coupling constant 
g, not even in the limit as a = fl4|]- It was demonstrated in |§ that, starting 



from a discretized, regularized version of the theory, the Euclidean theory (quan- 
tum Liouville theory) is obtained if the "bare" dimensionless coupling constant go 
is of order one. However, the relation between the bare coupling constant and the 
dimensionful continuum coupling constant g used in the present article is given 
by 

go = ga 3 , (7.1) 

where a is the lattice spacing in the dynamical triangulations providing the regu- 
larization (c.f. Sec. |2]). As discussed in | |14|| , the generalized CDT continuum limit 



corresponds to g fixed, a — > 0, and thus to go(a) — > 0. The fact that go(a) goes 
to zero in the CDT string field theory is of course related to the finite number of 
baby universes generated in this theory. By contrast, we have an infinite number 
of baby universes generated in non-critical string field theory, where go is of order 
one. 

However, there is clearly a deeper connection between the Euclidean and the 
CDT theory awaiting to be fully understood. It was shown in [|1^] that by inte- 
grating out the "excessive outgrowth" of baby universes in Euclidean 2d quantum 
gravity, one recovers the CDT theory, and the mapping between the dimensionless 
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variables x/V\ of the two theories was given explicitly. This mapping was later 
discovered by Seiberg and Shih p2 as the uniformization map from the algebraic 
surface representing the "semiclassical" non-critical string to the complex plane. 
The singular points of this algebraic surface correspond to so-called ZZ-branes, 



where there is a transition from compact to non-compact topology [23]. These 



singular points are mapped to points in the complex plane where one has a similar 
transition from compact to non-compact geometry in the CDT context |24 . 



It would be interesting to generalize the present string field theory based on 
causal dynamical triangulations to include the coupling to matter. In particular, 
one would like to investigate whether this theory still exhibits any trace of the 
presence of a c = 1 barrier. Since the existence of this barrier in the Euclidean 
theory can be partly understood as the result of an excessive creation of baby 
universes, tearing apart the two-dimensional worldsheet |5], 25], it is clear that 
the CDT theory may behave differently. Numerical simulations are compatible 
with the presence of a barrier at large values of the conformal charge c [26], but 
no definite results are available at this stage. Work is in progress on determin- 
ing whether the CDT string field theory can provide a useful analytic tool in 
addressing this situation. 

Equally interesting is the possibility of performing a summation over world 
sheets of all genera. Again, since the double-scaling limit in CDT string field 
theory is different from the double-scaling limit in non-critical string theory, and 
since there is a larger "penalty" for creating a higher-genus surface in the sense 
outlined above, viewing the creation of a higher-genus world sheet as a successive 
creation and annihilation of a baby universe, one could hope that the result of 
such a summation was better behaved and less ambiguous than was the case in 
non-critical string theory. Work in this direction is also in progress. 
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